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Abstract: We discuss the effective interactions of axion supermultiplet, which might be 
important for analyzing the cosmological aspect of supersymmetric axion model. Related 
to axino cosmology, it is stressed that three seemingly similar but basically different quan- 
tities, the Wilsonian axino-gluino-gluon coupling, the 1PI axino-gluino-gluon amplitude, 
and the PQ anomaly coefficient, should be carefully distinguished from each other for cor- 
rect analysis of the thermal production of axinos in the early Universe. It is then noticed 
that the 1PI axino-gluino-gluon amplitude at energy scale p in the range M$ < p < vpq 
is suppressed by M\jp 2 in addition to the well-known suppression by p/Wir 2 vpQ, where 
M$ is the mass of the heaviest PQ-charged and gauge-charged matter supermultiplet in 
the model, which can be well below the PQ scale vpq. As a result, axino production at 
temperature T > M$ is dominated by the production by matter supermultiplet, not by 
the production by gauge supermultiplet. Still the axino production rate is greatly reduced 
if M$ <C vpq, which would make the subsequent cosmology significantly altered. This 
would be most notable in the supersymmetric DFSZ model in which M$ corresponds to 
the Higgsino mass which is around the weak scale, however a similar reduction is possi- 
ble in the KSVZ model also. We evaluate the relic axino density for both the DFSZ and 
KSVZ models while including the axino production in the processes involving the heaviest 
PQ-charged and gauge-charged matter supermultiplet. 
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1. Introduction 

Supersymmetric extension of the standard model provides perhaps the most appealing 
solution to the gauge hierarchy problem [|J. Furthermore it can easily accommodate the 
axion solution to the strong CP problem ||, ||] which is another naturalness problem of the 
standard model. Supersymmetric axion model ||, || necessarily contains the superpartners 
of axion, the axino and saxion, which can have a variety of cosmological implications. 
In particular, axino can be a good candidate for cold dark matter Q in the Universe, 
depending upon the mechanism of supersymmetry breaking and the presumed cosmological 
scenario || [9| . 

One of the key issues in axino cosmology is the thermal production of axino by scatter- 
ing or decay of particles in thermal equilibrium in the early Universe. Most of the previous 



analysis of thermal axino production |T(], 11, 12] is based on the following form of local 



effective interaction between the axion supermultiplet and the gauge supermutiplet: 

d 2 6-^—W aa WZ, (1.1) 
32tt 2 v P q 

where vpq is the scale of spontaneous PQ breakdown, and A = (s + ia)/y2 + \[2Qa + 
6 2 F A is the axion superfield which contains the axion a, the saxion s, and the axino a 
as its component fields. If axinos were in thermal equilibrium, their relic density can be 
determined in a straightforward manner, and the result does not depend on the details 
of axino interactions. For the other case that axinos were not in thermal equilibrium, 
one usually assumes that the production is mostly due to the effective interaction (|1 . 1|) . 
Then a simple dimensional analysis implies that the axino production rate per unit volume 
scales as oc T 6 /(167t 2 wpq) 2 , which results in a relic axino number to entropy ratio 
Ya oc Tr/(16tt 2 vpq) 2 , where the reheat temperature Tr is assumed to be lower than both 
vpq and the axino decoupling temperature. If axino is a stable dark matter (DM), this 



- 1 - 



result leads to a severe upper limit on m^Tp ||, g|. Even when axino is not stable, thermal 
axino production in the early Universe can affect the later stage of cosmological evolution 
in various ways, e.g. there can be a late decay of axino into lighter particles in the minimal 



supersymmetric standard model (MSSM) or into the lighter gravitino and axion [13, 14 



15, 16], which might affect the relic DM density, or the Big-Bang nucleosynthesis, or the 



structure formation. 

In this paper, we discuss the effective interactions of axion supermultiplet in a general 
context, and examine the implications for cosmological axino production. It is stressed 
that three seemingly similar but basically different quantities, the Wilsonian axino-gaugino- 
gauge boson coupling, the 1PI axino-gaugino-gauge boson amplitude, and the PQ anomaly 
coefficient, should be carefully distinguished from each other for correct analysis of the 
cosmological axino production. At any given scale, the Wilsonian coupling is local by 
construction, but it depends severely on the field basis adopted to define the effective 
theory. On the other hand, the 1PI amplitude is an observable quantity, and therefore 
is basis-independent, while it generically contains non-local piece and has different value 
at different energy scale of external particles. The PQ anomaly coefficient is an intrinsic 
property of the PQ symmetry of the model, which is basis-independent and has a common 
value at all scales as a consequence of the anomaly matching condition. 

At any temperature below vpq, axino production by gauge supermultiplet is deter- 
mined by the 1PI axino-gaugino-gauge boson amplitude at the corresponding energy scale. 
Our key observation is that the 1PI axino-gaugino-gauge boson amplitude at energy scale 
p in the range M$ < p < vpq does not scale like the Wilsonian effective interaction (|1 . 1|) , 
but is suppressed further by the factor M|/p 2 , where M$ is the mass of the heaviest 
PQ-charged and gauge-charged matter field in the model, which can be well below vpq in 
general. This is a simple consequence of that the axion supermultiplet is decoupled from 
gauge supermultiplets (and from charged matter supermultiplets also) in the limit M$ — > 0, 
which is not manifest if one considers the effective interaction (1.1) alone, but recovered in 
the full analysis taking into account all interactions of the axion supermultiplet together. 
As a result, if M$ <C vpq, which is possible in most cases, axino production by gauge 
supermultiplets is greatly reduced. This would be most dramatic in the supersymmetric 
Dine-Fishler-Srednicki-Zhitnitsky (DFSZ) model [17] as the model does not contain any 
exotic PQ-charged and gauge-charged matter field other than the matter and Higgs fields 
in the MSSM, and therefore M$ is around the weak scale. A similar suppression is possible 
in the supersymmetric Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion model [^] also, 
in which M$ corresponds to the mass of an exotic PQ-charged quark multiplet, which in 
principle can take any value between the PQ scale and the weak scale. 

An immediate consequence of the above observation is that axino production in the 
temperature range M$ < T < vpq is dominated not by the production by gauge super- 
multiplets, but by the production by the heaviest PQ-charged and gauge-charged matter 
multiplet <£. If M$ <C upg, axino production rate at T > M$ is greatly reduced compared 
to the previous result based on the effective interaction (LI) alone. This can significantly 
alter the cosmological aspect of the model. For instance, a high reheat temperature which 
has been considered to produce too much axino dark matter can be allowed if M$ <C vpq. 
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Also, our result can relieve the upper bound on the reheat temperature, which has been 



obtained in [14] by considering the gravitino and axino productions together, and might 
have implication for the cosmological bounds on axion supermultiplet in the presence of 
small imparity breaking, which have been discussed in [16]. Since the interaction rate be- 



tween the axino and the thermal bath of gauge-charged particles is reduced in the limit 
M<j> <C T, the axino decoupling temperature can take a higher value than the previous 



estimate based on (1.1). Similar observations apply to the axion and saxion cosmology 



also, and our result might affect some of the results of the recent analysis of thermal axion 



production jl9|. 



The organization of this paper is as follows. In the next section, we discuss in a 
general context the effective interactions of axion supermultiplet which will be crucial for 
later discussion of cosmological axino production. In section 3, we provide an analysis 
of thermal axino production for both the KSVZ and DFSZ models, and determine the 
relic axino density including the contribution from the processes involving the heaviest 
PQ-charged and gauge-charged matter multiplet at T > M$. Section 4 is the conclusion. 

2. Effective interactions of axion supermultiplet 

In this section, we discuss the generic feature of the effective interactions of axion supermul- 



tiplet [20], which will be relevant for our later discussion of cosmological axino production. 
For the purpose of illustration, we first consider a simple specific model, a supersymmetric 
extension of the KSVZ axion model flql , and later generalize the discussion to generic 
supersymmetric axion models. 

At the fundamental scale M*, which is presumed to be of the order of the reduced 
Planck scale Mpi ~ 2.4 x 10 18 GeV, the supersymmetric part of our model is described by 
the following high scale lagrangian 



£(M*) = J d A 6K+ J d 2 9 (^f a W aa W% + W(^!) \ + h.c. 



(2.1) 



where W£ are the gluon superfields, and the Kahler potential K, the holomorphic gauge 
kinetic function f a , and the superpotential W are given by 

1 



K = J2^I, fa 



W = hZ (yXY - V -^\ + XXQQ C , (2.2) 

where gl{M*) is the holomorphic QCD coupling at M*, and stand for the chiral matter 
superfields in the model, = {Z, X,Y,Q,Q C }. Among these matter fields, Z, X and 

Y are gauge singlet, while Q + Q c are colored quark multiplet. Here we distinguish the 
holomorphic QCD coupling g s from the physical QCD coupling g s for later discussion. 
Also we ignore non-renormalizable operators suppressed by 1/M*, as well as the similarly 
suppressed supergravity effects, under the assumption that the PQ scale vpq is far below 
M*, so that any effect suppressed by vpq/M* can be safely ignored. 
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The model is invariant under the global PQ symmetry 1 

U(1) PQ : <&/ -> e ixia ^j, (2.3) 
where the PQ charge xi are given by 

(x z , x x , x Y , x Q + xqc) = (0, -1, 1, 1) . (2.4) 
This PQ symmetry is explicitly broken by the QCD anomaly as 

«Q = ^CpqF^F^, (2.5) 



where Jpn is the Noether current for the symmetry transformation ( |2.3D defined in a gauge 
invariant regularization scheme, and then the anomaly coefficient is given by 

C PQ = 2^x/Tr(T c 2 (<D 7 )) = 2(x Q + x Qc )Tr(T c 2 (Q)) = Nq, (2.6) 
I 

where T c (<f>/) is the generator of the SU(3) C transformation of and Nq is the number 
of Q + Q c . The QCD instantons then generate an axion potential which has a minimum 
at the axion vacuum expectation value (VEV) which cancels the QCD vacuum angle, so 
the strong CP problem is solved Q. In our convention, the PQ anomaly coefficient Cpq 
corresponds to the number of the discrete degenerate vacua of the axion potential, which 



can have a nontrivial cosmological implication [21]. In the following, we will assume Nq = 1 
for simplicity. 

In the limit to ignore SUSY breaking effects, the model has a super symmetric ground 
state 

(XY) = V -^, (Z)=0. (2.7) 

Including SUSY breaking effects, the VEV of the ratio X/Y can be fixed and a small but 
nonzero VEV of Z can be induced also. In fact, SUSY breaking part of the model is not 
essential for our discussion of axino effective interactions, although it is crucial for the 
determination of the axino and saxion masses [^, 23 [ . We thus simply assume that soft 
SUSY breaking terms in the model fix the VEVs as 

(X) ~ (Y) (2.8) 

with appropriate (model-dependent, but light) axino and saxion masses, for which the 
axion scale is given by vpq. Also, to illustrate our main point, we further assume that 



In addition to the PQ symmetry, we can introduce an approximate U(1)r symmetry under which 
8 a — ► e l ^0 a and Z — > e 21 ^ Z to justify the form of the superpotential. In the present form, the model is 
invariant under another t/(l) symmetry which makes the massive Dirac quark Q + Q c stable. However it is 
straightforward to break this (7(1), while keeping the PQ symmetry unbroken, when the model is combined 
with the MSSM, making Q + Q c decay fast enough. 
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h = 0(1), but A <C 1, so the model involves two widely separated mass scales, the axion 
scale vpq and the quark mass Mq which can be far below vpq: 



Mq = < vpq. 



(2.9) 



All physical consequences of the model can be determined in principle by the high scale 
lagrangian (2.1). However, for low energy dynamics at scales below vpq, it is convenient to 
construct an effective lagrangian at a cutoff scale A just below vpq, which can be obtained 
by integrating out the massive fields heavier than A as well as the high momentum modes 
of light fields at scales above A. In our case, two chiral superfields among X, Y and Z get 
a mass of 0(vpq), while the remained combination provides the axion supermultiplet. To 
integrate out the heavy fields, one can parameterize them as 

X = -L(v PQ + Pl )e- A ^, Y = -L(v PQ + Pl )e A / v ^, Z = p 2 , (2.10) 

for which the Kahler potential and superpotential at the UV scale M* are given by 

K = k; PQ + W | 2 cosh (^j^-) + P%Pi + Q ] Q + Q C] Q C , 

W = hv PQPlP2 + hp 2P j + (m q + V2^- Pl ) e~ A / Vp QQQ c , (2.11) 
2 V V PQ J 



where pi (i = 1, 2) are massive chiral matter fields, and 

A = -|= (s + ia) + V29a + 6 2 F A 
V2 



(2.12) 



is the axion superfield which contains the axion a, the saxion s and the axino a as its 
component fields. 

After integrating out the heavy p% and the high momentum modes of light fields, the 
resulting Wilsonian effective lagrangian can be chosen to take the following form 



£eff(A) 



K A (A + A*) + Z Q (A + At)Qtg + z Qc (A + A^Q C 
j d 2 9 (Iff {A)W aa W* + Wes) + h.c. 



(2.13) 



where 



K A = -(A + A^f + O 
1 C w A 



{A + At)- 

VPQ 



off 



fl 9lW tor 2 v PQ ' 
W eS = M Q e-^ +3: Q c)A/vpc iQQ c . 

With the above form of effective lagrangian, the PQ symmetry is realized as 

U(1) PQ : A^A + iavp Q , $ -> e^ a ^ ($ = Q,Q C ), 



(2.14) 



(2.15) 
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and the PQ anomaly in the effective theory is given by 



Cpq = C w + xq + xqc, (2.16) 



which should be same as the PQ anomaly (|2.6| ) in the underlying UV theory. For generic 
form of effective lagrangian consistent with the PQ symmetry ( 2.1 5| ) , there can be three 
types of axino effective interactions at the order of 1/vpq: 

A]£(A) = - / d 2 9^r—W aa WZ + h.c., (2.17) 

J J27T^ V P Q 

A 2 £(A) = / d 4 9 ^ A + A) (y Q QiQ + y QC Q^Qc) , (2.18) 

J VpQ V ' / 

A 3 £(A) = - [ d 2 6 (x Q + x c q )Mq QQ C + h.c, (2.19) 

J ' VpQ 



where 



91nZ$ 

y$ = vpq 



dA 

In our particular case, one easily finds 



($ = Q,Q C ). (2.20) 



A=0 




C w = 0, x$ = x*, m(A) = O —— 2~T~^ In T , (2-21 



where the small y$ is induced by the radiative corrections involving the Yukawa coupling 
A = \^2Mq/vpq at UV scales above A. This result shows that the axion supermultiplet 
is decoupled from the gauge multiplet and charged-matter multiplet in the limit Mq — > 0, 
which is manifest in the UV theory ( 2.11| ). Note that the PQ symmetry splits into two 



global U(l) symmetries in the limit Mq — > 0, the axial U(l) symmetry of Q,Q C and the 
anomaly free U(l) not involving the transformation of Q,Q C , and the axion becomes a 
massless Goldstone boson of the anomaly-free U(l) part in the limit Mq — > 0. 

In fact, one can consider a different description of the same effective theory, for instance 
a scheme using different form of PQ symmetry under which all light fields at A are invariant 
except the axion superfield. In our case, such description can be achieved by the A- 
dependent field redefinition <I> — > e x ® A / Vp Q&, making all matter fields except A invariant 
under the PQ symmetry. Such form of PQ symmetry can be convenient in certain respects 
since x$ = (<E> = Q, Q c ), so the PQ anomaly originates entirely from the variation of the 



local effective interaction (2.17), which results in Cpq = Cw- However such field basis is 
not convenient for a discussion of physics at energy scales above Mq since the decoupling 
of the axion supermultiplet from the gauge and charged-matter multiplets in the limit 
Mq — > is not manifest, but is achieved by the cancellation between the contributions 
from C w = N Q and yg + yQc = 1. 

To get more insights, let us consider a more general description associated with the 
field redefinition 

$ -> e cx * A/vp Q<$> ($ = Q, Q c ), (2.22) 
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where c is an arbitrary real constant parameterizing the field basis. After this field redefi- 
nition, the PQ symmetry is realized as 



U(1) PQ : A^A + iavp Q , 3> -> e^ 1 -*)**"*, (2.23) 
and there appears an axion-dependent local counter term in the effective lagrangian 

(2.24) 



/ 



32n z 



V PQ 



which is due to the Konishi anomaly [24| for the field redefinition ( 2. 22] ). Note that this 



Konishi anomaly term is required to match the PQ anomaly of (2.23) with the PQ anomaly 
J2.6| ). Including the Konishi anomaly, Wilsonian effective couplings of the axion superfield 
in the redefined field basis are given by 



C w = c(x Q + XQc), x$ = (1 - C)X0, 



y$(A) = cx$ + O 



16vr 2 v 2 pQ V A 



(2.25) 



Note that the PQ anomaly in the redefined field basis is determined by two contributions: 



CpQ = C(XQ + XQc) + (1 - c)(x Q + XQc), 



(2.26) 



where the first piece is from the variation of the local effective interaction ( p. 24 ), while 
the second piece is from the axial anomaly of the PQ transformation of Q,Q C in ( 2.23j ). 
Obviously each contribution depends on the field basis parameter c, but their sum is 
independent of c as it should be. 

Now, if one ignores the part of y$ suppressed by M^/vpq, Wilsonian effective inter- 
actions of the axion supermultiplet are given by 



Ai£(A) 
A 2 £(A) 
A 3 £(A) 



d*e c M±^- ^-w™wz,+h.c., 



32^ 



vpq 



(A + ^t) 
V PQ 



cxqQIQ + cxqcQ c ^Q c 



I d 2 6 (1 - c){x Q + x c Q )M Q —QQ c + h.c. 
J ' vpq 



(2.27) 
(2.28) 
(2.29) 



With these, one can compute the axino production rate at any temperature T < A, and 
the result should be independent of the field basis parameter c since all physical quantities 
should be basis- independent. However, most of the previous studies use only the effective 
interaction of the form ( 2.27| ), which can lead to a misleading result as we will see below. 
In fact, considering only the effective interaction (2.27) leads to a highly overestimated 
result for the axino production at energy scale p ^> Mq. At such high energy scale, one 
can take a limit Mq/p — > in which the axion supermultiplet is decoupled from the gauge 
multiplets and charged-matter multiplets. Although such decoupling is not manifest in the 
effective interactions ( |2,27| ) — ( 2.29[ ) with c 7^ 0, it is manifest in the field basis with c = 0, 
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(a) (b) (c) 

Figure 1: Diagrams for the process g + g —> g + a. 




(a) (b) (c) (d) 

Figure 2: Contributions to the 1PI axino-gluino-gluon amplitudes from the loops of Q, Q c . 



and also in the underlying UV theory ( 2.11 ). This means that axino production by gauge 
supermultiplets at p > Mq should be suppressed by some powers of Mq/p if one takes into 



account the interactions ( 2.27 ) — ( 2.29 ) altogether to get a correct c- independent result. 
On the other hand, the interaction ( p. 27 ) by itself does not involve any suppression by 
Mq/p, and therefore the analysis using ( pT27l ) alone gives a highly overestimated result in 
the limit p 3> Mq. 



More explicitly, the effective interaction (2.27) gives the Wilsonian axino-gluino-gluon 
amplitude 



A w (k,q,p) 



g 2 c(x Q + x Q c 
16n 2 y/2vpQ 



-5 4 (/c + q + p)u(k)a^v(q)e^p' y 



(2.30) 



where p^ and are the gluon momentum and polarization, and u{k) and v(q) are the 4- 
component Majorana spinor wavefunction of the axino and gluino, respectively. If one uses 
this amplitude alone to compute the axino production rate in the process gluon + gluon — > 
gluino + axino (see Fig. [l]), a simple dimensional analysis tells that the rate (per unit 
volume) is given by 



r(55 -> go) 



c 2 (xq + x Q c) 2 



£<7 S 6 T 6 



(l6ir 2 vpQ) 2 ' 



(2.31) 



where £ is a dimensionless coefficient which is independent of c. However this can not 
be the correct answer as it depends on the field-basis parameter c, and there should be 
additional contribution which removes the c-dependence of the result. Indeed, including 
the contributions from the axino-gluon-gluino amplitude due to the loops of Q, Q c (see Fig. 
|2|), which involve the c-dependent axino- matter couplings in ( 2.28| ) and ( 2.29j ), we find that 
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the c-dependence in the production rate disappears as required, and the final result takes 
the form 

r(gg ga) = i 2 rM^, ( 2 - 32 ) 



(1§-k 2 v P q) 2 ' 



where 



^°(#'" 2 (i^)) lbr T>M '- 



7 = {xq + 4) + O for T < M Q , (2.33) 

which shows that the production rate is indeed suppressed by Mq/T 2 at T > Mq. This 
result can be obtained by replacing the axino-gluino-gluon vertex in Fig. [I] with the one- 
particle-irreducible (1PI) axino-gluino-gluon amplitude 

A 1PI (k,q,p) = - P) 5\k + q+ p)u(k)W5v(q)e»p» (2.34) 

which has the following limiting behavior (see the appendix A) when the axino and gluino 
are on mass-shell: 



C 1PI {k 2 = q 2 = 0; \p 2 \ » Mq) = (xq + x Q )^§ [in 2 j + O In 



\P 2 \ 



ClPI (k 2 = q 2 = 0; \p 2 \ « Mq) = Xq + Xq + ^ j . (2.35) 

Generically the 1PI amplitude due to light particle loops includes non-local and non- 
analytic piece in the limit p 2 3> M 2 ht , while the contribution from heavy particle loop 



allows a local expansion in the limit p 2 <C M? . For the 1PI axino-gaugino-gauge boson 



'light' 

heavy ' 

amplitude, it takes the following form (see the appendix A) when the axino and gaugino 
(or gauge boson) are on mass-shell: 



where Cipj can have a logarithmic p-dependence due to higher loop effects. The constant 
(or logarithmic) C\ P i can be encoded in the PQ-invariant real 1PI gauge kinetic function 
J- a which defines the 1PI gauge kinetic term 26]: 

[ d 4 e^-T a (lnp 2 ,A + A^)W aa ^W^ (2.37) 
J 16 p z 



and then Cf PI is determined as 



2 



C 1P i = -8tt v pq — 



(2.38) 

A=0 
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This 1PI operator should be distinguished from the Wilsonian gauge kinetic term 



/ 



d 2 0jff (In A, A)W aa W% + h.c. (2.39) 



which determines Cw as 



a., dff 



C w = -8vr v P Q 



(2.40) 

A=0 



As we have noticed, Cw is a basis-dependent Wilsonian coupling which changes under an 
yl-dependent holomorphic redefinition of charged matter fields, while C\pi is an observable 
amplitude invariant under the redefinition of matter fields. 

In the above, we find in the context of a simple supersymmetric KSVZ axion model that 
the 1PI axino-gaugino-gauge boson amplitude at energy scale in the range M$ < p < vpq 
is suppressed by M\jp 2 in addition to the well-known suppression by p/16tt 2 vpq, where 
M<j> is the supersymmetric mass of the heaviest PQ-charged and gauge-charged matter 
field, which can be well below vpq. On the other hand, the 1PI amplitude at energy scales 
below M$ has the same scaling behavior as the Wilsonian amplitude ( |2.30 ), which can 



be noticed from (2.35), although the precise coefficient can be different in general. This 



observation applies also to another type of popular axion model, the DFSZ model [17|. To 
see this, let us consider a supersymmetric extension of the DFSZ axion model described 
by the following form of the Kahler potential, gauge kinetic function and superpotential at 
the UV scale M*: 

I 



( v 2 \ X 2 

W = hziXY-^\+K jr H u H d + --- , (2.41) 

where denote the matter fields in the model, including the MSSM Higgs, quark and 
lepton superfields, and the ellipsis stands for the PQ-invariant Yukawa couplings for the 
quark and lepton masses. The matter fields transform under the PQ symmetry as 

U(1) PQ : -> e* 8 ' '*/ (2.42) 

with the PQ charges 

(xz,x x ,XY,x Hu ,x Hd ) = (0,-1,1,1,1), (2.43) 

and the PQ charges of the quark and lepton superfields can be fixed by the MSSM Yukawa 
couplings. Here the Higgs /x-term is generated by the spontaneous breakdown of the PQ 
symmetry || |7|, yielding 

2 

_ ^PR_ ( 2 44) 

The most notable feature of this DFSZ model is that the Higgs doublets correspond to the 
heaviest PQ-charged and gauge-charged matter field, so M<j, = [i which should be around 
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the weak scale. It is straightforward to repeat our discussion for this DFSZ model, and 
then one finds the 1PI axino-gaugino-gauge boson amplitude at p > \i is suppressed by 
fj, 2 /p 2 . This would lead to a dramatic reduction of the cosmological axino production, when 
compared to the result of the previous analysis using only the effective interaction of the 
form ( |2.17 ). It should be noted that a similar reduction is possible for the KSVZ model 



also as M$ = Mq can be in principle comparable to the weak scale. 

Our result can have important implications for the cosmology of supersymmetric ax- 
ion model. For instance, certain parameter range of the model and/or the cosmological 
scenario, which has been considered to give a too large relic axino mass density in the pre- 
vious analysis, can be safe if the model is assumed to have M$ <C vpq. It also implies that 
at temperature T > M$, axinos are produced dominantly by the processes involving the 
heaviest PQ-charged and gauge-charged matter supermultiplet, e.g. Q + gluon — y Q + a for 
the KSVZ model and higgs + W — > higgsino + a for the DFSZ model, since the amplitude 
of such process does not involve the loop suppression factor and is suppressed only by a 
single power of M^/vpq. In the next section, we compute the axino production rate and 
the resulting relic axino number density while including the production by the heaviest 
PQ-charged and gauge-charged matter multiplet for both the KSVZ and DFSZ models. 

It is in fact straightforward to generalize our discussion to generic supersymmetric 
axion models. For this, let us consider a general form of Wilsonian effective lagrangian of 
the axion superfield at a scale A just below vpq, which takes the form 



C cS (A) = J d A 9 (K A (A + rf) + Z n {A + tf)& n $ 
f d 2 9 (^ff(A)W aa W« + W cS ^ +h.c 



(2.45) 



where {& n } denote the light gauge-charged matter fields at A, and 
2 \ vpq J 



V PQ V Vp Q 



./ a 



Hi _ 1 @W A 



02(A) 8^ v PQ ' 



W eS = l -e-^ + ^ A ' Vp <iM mn <5> m <5> n + ^e-^ + ^ +5: ^ A l Vp Q\ mnp <$> m <5>n%. (2.46) 
The PQ symmetry in this effective theory is realized as 

U{l) PQ : A^A + iavp Q , $ n -> e^ Q $ n , (2.47) 
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and the Wilsonian effective interactions of the axion superfield are given by 

C« r A 



Ai£(A) = - / d 2 6 — %r W aa W® , (2.48) 

J 32tt z v P q 

A 2 £(A) = f d A 9 y n ( A + A) *t* n , (2.49) 
J V PQ 

/A r 1 
d 2 6 -(x m + x n )M mn <$> n <S> m 

vpq L 2 

+ -(z m + x n + x p )A mnp $ m $ n $p . (2.50) 

According to our discussion, there are three quantities which are all related to the axino 
coupling to gauge supermultiplets, but basically different from each other: 

{ C{y> CpQ, Cipj }, (2-51) 



where are the Wilsonian couplings in ( 2.4E ), C P q are the PQ anomaly coefficients 
defined as 



«q = ^C a PQ F^F^, (2.52) 



and finally Cf PI determines the leading part of the 1PI axino-gaugino-gauge boson ampli- 
tude 



,2 



A a lPI (k,q,p) = -—f- CZ PI 6*(k + q+p)u(k)°^v(q)e^ (2.53) 

167T z V2fpQ 

which shows the behavior 

C? PI (k 2 = q 2 =0, Mi ght < p 2 < M h 2 eavy ) 

where Mught and Mh eav y denote the masses of matter fields in the effective theory ( 2.45| ). 
It is then straightforward to find 

C PQ = C^ + 2^5 n Tr(T a 2 (cI> n )), (2.55) 

n 

C a lPI {M 2 <p 2 < A 2 ) = Cfr - 2^y n Tr(T a 2 (<D n )), (2.56) 

n 

where M$ is the mass of the heaviest PQ-charged and gauge-charged matter field in the 
model. Here the expression of C P q is exact and valid at any scale below vpq, while the 
expression of C^ PI is derived in 1-loop approximation. It is also straightforward (see the 
Appendix A) to compute Cf PI at lower momentum scale, which yields 

C a 1PI {p 2 <M 2 ) =Cfr(A)-2 Y, yn(A)Tr(T 2 (cD n )) 

+ 2 Y 5m(A)Tr(T a 2 ($ m )) (2.57) 

Ml>p 2 
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in 1-loop approximation. 

In fact, one can easily derive an exact (in perturbation theory) expression of C^ PI , 
including the piece depending logarithmically on the momentum scale p. We already 
noticed that Cf PI can be determined by the real 1PI gauge kinetic function F a as 



Cipj = -8tt 2 vpq 



OA 



(2.58) 



A=0 



Solving the Novikov-Shifman-Vainshtein-Zakharov RG equation in the limit M|/p 2 — > 0, 

Tr(T 2 (G<)) 



one finds 25, 26] 



T a = Re(ff) + 



16vr 2 



hi 



El 

A 2 



8tt 2 



8tt 2 



lnJ- a , (2.59) 



where b a = — 3Tr(T 2 (G)) + ^ n Tr(T a 2 ($„)) is the one-loop beta function coefficient, and 
7 n = dlnZ n /dlnp 2 is the anomalous dimension of <£ n for the 1PI wavefunction coefficient 
Z n which can be chosen to satisfy the matching condition 



Z n (p 2 = A 2 ) = Z n (A). 



We then find 



where 



Ct PI {Mi < p l < A 2 ) 



l-Tr(T 2 (G)) 5 2 (p)/8^ 



Vnip) = V P Q 



dlnZr, 



dA 



(2.60) 



(2.61) 



(2.62) 



A=0 



which gives (2.56) at 1-loop approximation. 

Within the effective theory ( 2,45| ), one can make a holomorphic redefinition of matter 
fields 



after which the PQ symmetry takes the form 

U(1) PQ : A -> A + iavpg, $ n -> e^"" 2 "^, 
and the Wilsonian couplings of the axion superfield are changed as 

C^C^ + 2j> n Tr(T 2 ($ n )), 



(2.63) 



(2.64) 



Vn 



(2.65) 



This shows that one can always choose a field basis with x n = 0, for which only the 
axion superfield transforms under the PQ symmetry, and therefore = C P q. Another 
interesting choice would be the field basis with y n = 0, for which Cyy = C^ PI . Note 
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that Cpq and Cf PI are directly linked to observables, and therefore invariant under the 
reparametrization ( |2,65 ) of the Wilsonian couplings associated with the field redefinition 



(2.63). It should be noticed also that for a given theory, C P q have common values at all 
scales, while Cf PI can have different values at different momentum scales. On the other 
hand, are lagrangian parameters which can take different values in different field basis 
(or in different UV regularization) within the same theory. 

A key result of our discussion, which has direct implication for cosmological axino 
production, is that the 1PI axino-gaugino-gauge boson amplitude in the momentum range 
M§ < p 2 < vpg is suppressed by M|/p 2 , more specifically 2 

M 2 / M 2 



( M 2 

C a lPI {M 2 <p 2 < v P q) = C a 1PI {M 2 <p 2 < v 2 PQ ) + 0[-^-hx 



P 2 



M 2 , o (M 2 



(2 ' 66) 

where M$ is the supersymmetric mass of the heaviest PQ-charged and gauge-charged 
matter field in the model, which can be well below vpq in most cases. As we will see, there 
can be contributions to the above 1PI amplitude from UV dynamics at scales above vpq, 
but they are generically of 0{M% \n(M 1r /v PQ )/S'n 2 v 2 PQ ) or 0(v 2 PQ /M 2 ), where M* is the 
fundamental scale of the model, e.g. the Planck scale or the GUT scale, which is presumed 
to be far above the PQ scale. These UV contributions are either smaller than ( 2.66| ) or 



negligible by itself when M* is comparable to the Planck or GUT scale. As we will argue 
below, the result ( p. 66 ) applies to generic supersymmetric axion model if the model has a 



UV realization at M*, in which (i) the PQ symmetry is linearly realized in the standard 
manner, 

U(1) PQ : $ 7 -> e*** a * J} (2.67) 

where stand for generic chiral matter superfields, and (ii) all higher dimensional 

operators of the model are suppressed by appropriate powers of 1/M*. 

To proceed, let {&a} denote the gauge-singlet but generically PQ-charged matter fields, 
whose VEVs break U(l) pq spontaneously, and {$ ra } denote the gauge-charged matter fields 
in the model. Then the Kahler potential and superpotential at the UV scale M* can be 
expanded in powers of the gauge-charged matter fields as follows 

K = K PQ (^ A , * A ) + (l + ^t^s + • • • ) *t + • • • , 
W = Wpq^a) + l(\A mn $A + ^p^$A$B + • • -)*m*», 

+ \ (Xmnp + ^P$A + • • • ) + • • • , (2.68) 



2 This is for the kinematic regime with a gauge boson (or gaugino) 4- momentum having \p 2 \ > M%, 
while the axino and gaugino (or gauge boson) 4-momenta are on mass-shell. Cf P/ can have a bit different 
behavior in other kinematic regimes, but always suppressed by M|,/p 2 if any of the external particles has a 
4-momentum p with \p 2 \ > M%. For instance, if any of the external particles has a vanishing 4- momentum, 
while the other particles have |p 2 | > AfJ, the amplitude is of O \ -jjt m (p 2 /-Ml! 
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where Kpq and Wpq are the Kahler potential and superpotential of the PQ sector fields 
{3>a}j an d the ellipses stand for higher dimensional terms. Here we assume that the Kahler 
metric of the gauge-charged matter fields <£ n is flavor-diagonal for simplicity. Also, to be 
complete, we includes the leading higher dimensional operators suppressed by 1 /M* . Under 
the assumption that Kpq and Wpq provide a proper dynamics to break the PQ symmetry 
spontaneously, we can parameterize the PQ sector fields as follows 



(2.69) 



where va = (&a) with v P q = Y1a x a\ v a\ 2 , Pi denote the massive chiral superfields in the 
PQ sector, and U a% are the mixing coefficients which are generically of order unity. For 
this parametrization of the PQ sector fields, the Kahler potential and superpotential at 
M* take the form 

K = K PQ ( P l Pi , A + Ai) + (zg> +0(^(A + At), *t* n + , 

W = W PQ ( Pi ) + \ (M mn + e-^ + ^ A ^Q^ n 

+ \ (\ mnp + °(j£j) e-^ +x ^ A ' v ^^ n % + • • • , (2.70) 
where and M mn are field-independent constants, and the Yukawa coupling constants 



y mnp 



obeys the PQ selection rule 



{%rn %p) ^mnp — (^m ~i~ %n ~i~ -Bp) ^^mnp ~i~ Q f ^ j J — Q ( ^ ^ • (2-71) 

One can now integrate out the massive pi as well as the high momentum modes of light 
fields, and also make the field redefinition ( |2.63| ) to derive an effective theory in generic 
field basis. The resulting effective lagrangian at the scale A just below vpq takes the form 
of (ggp with 

Cfr = -S^vpq-Ll. = 2^, n Tr(T a 2 ($ n )), 

n 

- m _„ dlnZ n ( (x m + x n ) Ml n [MA\ , n (vp Q \ 

+ °(w)< (2 - 72) 

where we have set 1ii(M*/upq) = 0{ir 2 ) and used {x m +x n )M mn < 0(M$) for M$ denoting 
the mass of the heaviest PQ-charged and gauge-charged matter field in the model. Here 
Cyy and the first part of y n are due to the field redefinition (|2.63| ), the second part of 
y n is from the loops involving the Yukawa couplings between the PQ sector fields and 
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the gauge-charged matter fields, which are generically of 0(M mn /vpQ) and depend on the 
axion superfield through the combination (x m +x n )A, and finally the last part of y n is from 
the higher dimensional operator in the Kahler potential of There can be additional 
contribution to y n from the loops involving the Yukawa couplings X mnp obeying the PQ 
selection rule ( p. 71 ), which is still within the estimate of y n in ( |2.72| ). We then have 



C a 1PI (p 2 = A 2 ) = C^(A)-2^y n (A)Tr(T a 2 (d> n )) 



V PQ, 



at the cutoff scale A just below v P q, and the PQ selection rule ( 2.71 ) takes the form 

% 

J PQ, 



(x m +x n + x p + y m + y n + y p )\mn P = °[ ~JT~ j + ° (j^^j ' {2.74) 



In the Appendix B, we examine the IPI RG evolution of C® PI including higher loop 
effects, and show that the above estimate of Cf PI is valid at generic momentum scale in 
the range M$ < p < v P q. This implies that C^ PI in the momentum range M$ < p < v P q 
is indeed dominated by the piece of O (^-ln 2 (p 2 /M 2 )), so the estimate of Cf PI 

is valid even when higher loop effects are taken into account. We thus conclude that in 
generic supersymmetric axion model with a PQ scale hierarchically lower than the UV 
scale M*, which is presumed to be around the Planck scale or the GUT scale, the IPI 
axino-gaugino-gauge boson amplitude at momentum scales in the range M$ < p < vpq 
is suppressed by M|/p 2 , in addition to the suppression by p/16ir 2 v P Q, where M<j> is the 
mass of the heaviest PQ-charged and gauge-charged matter field in the model. With the 
boundary condition ( 2.73| ) at p 2 = A 2 , one can determine C± PI at lower momentum scale 



p < M$ by computing the threshold correction. Using the result obtained in the Appendix 
A, we find the leading constant part of C^ PI at generic momentum scale below v P q is given 
by 

C? PI ( P ) = C7^(A)-2 yn(A)Tr(T 2 (d> n )) + 2 £ £ m (A)Tr(T a 2 ($ m )) 

= 2 (x m (A)+y m (A))Tr(T 2 (<D m )), (2.75) 

M2 l >p 2 

where Cy V (A),y n (A) and x n (A) are the Wilsonian couplings in the effective lagrangian 
(2.45) at the cutoff scale A just below v P q. Note that this general result correctly repro- 



duces the IPI amplitude ( 2.35| ) at p < Mq in the KSVZ model. 



3. Thermal production of axino 

In this section, we examine the thermal production of axino with the effective interactions 
which generically take the form (|2.48|) — (|2.50 ). As we have noticed, if the model has a UV 



completion (at a scale M* 3> v P q) in which the PQ symmetry is linearly realized in the 
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standard manner and all non-renormalizable interactions are suppressed by the powers of 
1/M*, the effective interactions (2.48) — ( |2.50| ) are constrained by the matching condition 



d2.73|) at the scale A just below vpq. Then one can choose a field basis in which the 
Wilsonian couplings of axion supermultiplet at A are given by 

CW(A) = 0, x n = x n , y n (A) = O O^t, V -^\ . (3.1) 

Of course, one can choose any other field basis, for instance the one with 

C w = 2^2 ^Tr(T a 2 ($ n )), ~x n = x n - z n , y n = z n + 0\ -gf , (3.2) 



V PQ 



which would be obtained by the field redefinition ( 2.63 ). Then the above Wilsonian cou- 
plings for arbitrary real values of {z n } should give the same physical results as those of 
( |3.1| ). The field basis (3.1) is convenient for describing the physics at energy scales above 



M$ since the decoupling of the axion supermultiplet in the limit M$ — > is manifest. How- 
ever, for physics at lower energy scales below M$, it is often more convenient to choose the 



field basis (3^) with x n = x n — z n = for which Cw = Cpq. 

Let 3>,<I> C denote the heaviest PQ-charged and gauge-charged matter superfield with 
a supersymmetric mass M$. In the field basis (|3.1|), the relevant effective interaction of 
axion supermultiplet takes a simple form 

- I d 2 8(x$ + x< s> c)M$—<S»S> c + h.c., (3.3) 

J V P Q 

where we have ignored the small y n = 0{Ml/vp Q ,v PQ /M 2 ). In the KSVZ model (|J), 
<!>, <I )C correspond to an exotic vector- like quark multiplet with x$ + x$c = 1, while in the 



DFSZ model ( 2.41 ), <E>, <1> C correspond to the Higgs doublet superfields H u , Hd with M$ = /j, 
and x$ + x$c = 2. A key element for the axino production by gauge supermultiplet is the 
1PI axino-gaugino-gauge boson amplitude which is given by 

A lP i{k,q,p) = - 16?r2 ^ S4 ( k + g + P)CiPi(k, q,p)u{k)a av ^v{q)e^{p)p\ (3.4) 



with 



ClP1 (k 2 = q 2 = 0; p 2 » M 2 ) ~ (z« + x$c)-f In 2 ^ ) (3.5) 

CiPi(k 2 = q 2 = 0; p 2 « Ml) = x ^ +x ^ + 0\J-^j. (3.6) 

With the above 1PI amplitude and also the axino-matter coupling ( |3.3| ), we can calcu- 
late the thermal production of axinos in the temperature range of our interest. Following 
1 11], here we consider the axino production processes listed in Table |]. (See Fig. § — @ 



for corresponding Feynman diagrams.) Among these processes, the processes A, B and F 
produce axino through the 1-loop transition g — >■ g + a (or 5—7- g + a) whose amplitude 
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is given by the 1PI amplitude (3.4). On the other hand, other processes produce axino 
through both the tree-level transition <I> — > $ + a (or + a) and the 1-loop transition 

g — > g + a (or g — > g + a). To compute the amplitudes of these axino production processes, 
we will use the field basis ( |3.1[ ) for which the decoupling of the axino in the limit M$ -4 
is manifest. 

The 1PI amplitude ( |3.5| ) implies that the amplitude of the axino production through the 
transition g — > g + a in the temperature range M$ <C T < vpq is suppressed by M|/T 2 . As 
a result, in this temperature range, axinos are produced mostly by the transition <I> — > $ + a 
(or + a) with an amplitude cx (x$ + x$c)M$/vpq, and then the production 

rate is given by 

r a (M$ «T< v PQ ) = 5^(<r(J + J -4 o + • • • Mn/nj 

= 0(1) X (x 4 + x $ c) 2 y f , (3.7) 

where nj is the number density of the I-th particle species in thermal equilibrium. On the 
other hand, at lower temperature T <C M$, the matter multiplet <!> is not available anymore, 
and axinos are produced either by g — > g + a or by q — > q + a, where q denotes a generic 
light matter multiplet with M q < T. For the temperature range 8ir 2 M q < T -C M$, 
looking at the magnitudes of the involved transition amplitudes, one easily finds that 
axinos are produced mostly by the transition g — > g + a (or g — > g + a) with an amplitude 
■Agga (z* + x^,c)/1Qtt 2 vpq, which results in 

T~ a {8n 2 M q < T « AT*) = 0(1) x (a* + x^f fj\ , (3.8) 



where the 0(1) factor includes the thermal field theoretic effects discussed in |T! 

Solving the Boltzmann equation, the relic axino number density over the entropy 
density can be determined as (see the Appendix C) 

n a (T ) _ [ T « dT T- a 
Ya[l0) - s(T ) ~ J To T s(T)H(T) 

where Tp is the reheat temperature, s(T) = 2i: 2 g*T^ /45 is the entropy density, and H(T) = 
y/ Tt 2 g*/90T 2 /Mpi is the Hubble parameter for the effective degrees of freedom g* and the 
reduced Planck mass Mp\ = 2.4 x 10 18 GeV. We then find 

Ya(87T 2 M q <T R <^ Af*) = 0(1) x (a* + x^f ff^ 1 T R , (3.10) 
Y~ a {M* <.T R <. vpq) = 0(1) x (a* + x^) 2 " A M 2 Pl M 9 , (3.11) 

ZTT VpQ 



where g = 135\/T0/27r 3 g^ / ' 2 . Here we used the result of [12] for the first result, while the 
second result is derived in the Appendix C. 
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Process 



Feynman diagrams 



\M\ 2 {87r 2 M q <^T <M,j>) 



\M\ 2 (M^ ^T^vpq) 



A 
B 

C 
D 

E 
F 
G 
H 
I 
J 



g + g-> a + g 
g + g^ a+g 
q + g -+ a + q 

q + g^- a + q 
q + q-f a + g 
g + g^- a + g 
q + g^- a + q 
q + g -+ a + q 
q+q^a+g 
q + q^-a + g 



Fig. | 

crossing of A 
Fig. | 

crossing of C 
crossing of C 
Fig. | 
Fig. | 

Fig. 
crossing of G 
crossing of H 



4Ci(s + 2t + 2t 2 /s) 
-4C 1 (t + 2s + 2s 2 /t) 

2sC 2 
-2tC 2 
-2tC 2 
-8d{s 2 +t 2 + u 2 ) 2 /stu 

-AC 2 (s + s 2 /t) 
-2C 2 {t + 2s + 2s 2 /t) 

-AC 2 (t + t 2 /s) 
2C 2 (s + 2t + 2t 2 /s) 



C(4 
C(2- 
C(4 
C(2- 



suppressed 
suppressed 

suppressed 

2Afj . 2A/| 
+ 1 — 771 -\ ~ 



-M% 1 TmJ) 
3M| s-3Mg - 



+ 



2Ml 



t 



--Ml + TMl> 
t-3M| _ u-3Ml \ 



2M| 



Table 1 

.,61 fabc|2 



Processes of axino production. Here C 



%g 2 Ml\T l3 {<S>y\ 2 /v 2 



Ci 



g b \f abc \ 2 /256TT 4 v 2 pQ , and C 2 = 5 b |Xg (q)| 2 /2567r 4 w|, Q , where $ denotes the heaviest PQ-charged 
and gauge-charged matter multiplet with a gauge-charge matrix given by T™ ($). For T <C M$, q 
stands for generic gauge-charged matter with M q < T, while it means $ for M$ -C T <C wpq- 




(a) 



(b) 



(c) 



(d) 



Figure 3: Diagrams for the process A. Here black dots represent the 1PI axino-gluino-gluon 
amplitude. 




Fig. |8] summarizes our discussion of the relic axino density. It shows Y% oc Tr for 
Tr < 0.1M$, which is due to that in the temperature range &-K 2 M q < T < Q.1M&, axinos 
are produced mostly through the transition g -+ g + a (or g -+ g + a) with an amplitude 
■Agga = 0(T/1Qtt 2 vpq). (Here M q corresponds to the mass of the next heaviest PQ-charged 
and gauge-charged matter multiplet.) If one uses only the effective interaction of the form 
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(a) (b) (c) 

Figure 5: Diagrams for the process F. 




(a) (b) (c) 

Figure 6: Diagrams for the process G. 




Figure 7: Diagrams for the process H. 

32 j- / d 2 9AW a W a to evaluate the axino production by g — > g + a (or g — > g + a), as one 
did in most of the previous analysis, one would get oc Tr even for higher Tr > 0.1M$, 
as represented by the dashed line. Taking it into account that the 1PI axino-gluino-gluon 
amplitude A g¥a = 0{M% ln 2 (T/Mq>)/lQ-K 2 Tvp Q ) 3 for T > M$, and therefore the axino 
production at T > M$ is mostly due to the transition <1> — > + a (or <I> — > <I> + a) 
with an amplitude = 0(M$/vpq), one can easily understand the behavior of 

for Tr > 10M$, which is nearly independent of Tr. Note that the dashed line crosses 
the correct solid line at Tr ~ 10 3 M,j,, implying that the previous analysis based on the 

effective interaction f d 2 9AW a W a alone gives rise to an overestimated axion relic 

density for the reheat temperature Tr > 10 3 M$, while it gives an underestimated for 

0. 1M* < Tr < 10 3 M $ . 

We now proceed to apply our results to the two well-known types of axion models, 

1. e. KSVZ |I| and DFSZ @ models. In DFSZ model, the heaviest PQ-charged and 

3 This estimate is based on the simple identification p 2 ~ T 2 , which might not work for instance for 
the t-channel processes with the final state axino momentum parallel to the initial state gluino (or gluon) 
momentum |28). However, including the thermal mass of intermediate state gluon (or gluino), which is of 
O(gT), this simple approach does work for the purpose of the order of magnitude estimation. 
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Figure 8: Relic axino number density over the entropy density vs the reheating temperature 
T R (solid line). The dashed line represents the result one would get by using only the effective 
interaction of the form ^-4 f d 2 9AW a W a . 

gauge-charged matter can be identified as the MSSM Higgs doublets, and then M$ = fj, is 
around the weak scale. On the other hand, in KSVZ model, <£, <£ c correspond to an exotic 
quark multiplet with M$ which can take any value between the weak scale and vpq. Since 
our results can have more important cosmological implication when M^/vpq is smaller, 
we first discuss the case of DFSZ model. 

3.1 DFSZ axion model 

To apply our results to the DFSZ model in which <£, <£ c correspond to the MSSM Higgs dou- 
blets H u ,Hd, we choose M$ = /i = 10 3 GeV, and consider the axino production processes 
listed in Table |l| while identifying (g,g) and (<&,<&) as the SU(2)i gauge supermultiplet 
and the Higgs supermultiplets, respectively. We use g 2 = g\(J = 10 4 GeV) ~ 0.5, and the 
effective degrees of freedom g*(T > fj,) = 228.75. The results are depicted in Fig. |9[ In the 
left panel, we plot Ya as a function of the reheat temperature Tr for three different PQ 
scales, vpq = 10 11 GeV (solid), 10 12 GeV (dashed), and 10 13 GeV (dotted). Our result 
shows that has approximately a constant value for Tp > 10 4 GeV, which is given by 

Y & (T R > 10 4 GeV) ~ 7.6 x 1(T 8 ( 1Ql2GeV V ( 3 . 12 ) 

V V PQ J 

In the figure, the black dashed line represents the axino number density when axinos were 
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Figure 9: (Left) Axino number density vs reheat temperature Tr for vpq — 10 11 GeV (solid), 
10 12 GeV (dashed), and 10 13 GeV (dotted). (Right) Contours giving fl~ a h 2 = 0.1. 



in equilibrium with the thermal bath of gauge-charged particles, which is given by 

Y~ eq = — ~ 1.8 x 10- 3 . (3.13) 

9* 

As we can see from Fig. ^, the axino interactions with gauge-charged particles in DFSZ 
model are too weak for axinos to be in thermal equilibrium even for high reheat temperature 
comparable to vpq. 

In the right panel of Fig. |9[ we plot the contours giving the relic axino abundance which 
would explain the observed DM density Qdm^ 2 — 0.1 J?]]. As an approximate guideline, we 
can refer to hot, warm, and cold axino dark matter for axino mass in the range uij < 1 keV, 



1 keV < rria < 100 keV, and m„ > 100 keV, respectively [11]. Then as shown in the figure, if 
vpq > few x 10 11 GeV, axino can provide the cold dark matter in our Universe with correct 
relic abundance. In our case, there is no upper bound on the reheat temperature coming 
from the relic dark matter density, and therefore we can avoid the previous conclusion that 
cold axino dark matter scenario is viable only when the reheat temperature is relatively 



low as Tr < 10 GeV [11]. In this analysis, we have ignored the effects of SUSY breaking 
and electroweak symmetry breaking on the axino effective interactions, which should be a 
good approximation for T S> 10 3 GeV. We note that the Higgsino (or the Higgs boson) 
decay into axino can be sizable if one includes SUSY breaking effect around the weak scale. 
If one includes the electroweak symmetry breaking effect, stop decay also can give a sizable 



contribution [15]. At any rate, as long as Tr > M$, the axino production by the decay 
processes does not significantly alter our result depicted in Fig. [j| 

For axinos to be a successful cold dark matter, we should also consider the constraints 
from the big bang nucleosynthesis (BBN). The long-lived next lightest SUSY particle 
(NLSP) such as the bino-like or wino-like lightest neutralino x or the lighter stau f can be 
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T R in GeV m a in GeV 

Figure 10: (Left) Axino number density vs reheat temperature Tr for vpq = 10 11 GeV (solid), 10 12 
GeV (dashed), and 10 13 GeV (dotted) with M Q /v PQ = 1CT 5 . (Right) Contours giving n d h 2 = 0.1. 



problematic in BBN if it destroys the light primordial elements by emitting an energetic 
photon or hadronic shower through the decay into axino. Neutralino decay to gauge boson 
and axino, \ ~ ► l/% + 5 can be induced by the 1PI axino-gaugino-gauge boson amplitude 
( |3.4| ), and its rate is estimated as 

1 ( d V M 1 1 / M% \ 3 / 10 12 GeV y 



where gi is the SU{2)l coupling constant. On the other hand, the stau decay f — > r + a 
is induced by tree- level process in DFSZ model since there is axino-gaugino mixing due to 
the electroweak symmetry breaking, and its decay rate can be estimated as [29] 



™? #2<Lk ml „„ An 1 f m f \/10 12 GeV\ 2 

r(f -^T + a) — U2 wcak cos 2 9 W s n 4 3 =— ( ^ — 

V ' 16vr v P Q M 9 2 W H 10" 5 sV200 GeVy V V PQ J 



(3.15) 

Hence, the NLSP neutralino lifetime is 0(1(T 7 ) - ©(lCT 1 ) s, while the NLSP stau lifetime 
is ©(lCT 11 ) - 0(1(T 5 ) s for v PQ = 10 9 - 10 12 GeV. In order to avoid the BBN constraints, 
the lifetime of such long-lived particles is required to be shorter than 10 2 s [30], which can 



be easily satisfied in DFSZ model over a reasonable parameter range of the model. 
3.2 KSVZ axion model 

Let us now consider the KSVZ axion model in which the heaviest PQ-charged and gauge- 
charged field is an exotic quark multiplet Q, Q c with Mq which can take any value between 
the weak scale and the PQ scale. Fig. || suggests that the previous analysis of axino produc- 
tion based on the effective interaction 3 27r a„ / d 2 9AW a W a alone can be applied only for 
Tr ~ O.IMq. More specifically, axinos are produced mostly by the gluon supermultiplet at 
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T < O.IMq, while at higher temperature axino production is mostly due to the transition 
Q — > Q + a (or Q — > Q + a). Thus our discussion in this paper can significantly alter the 
previous result if Mq is well below vpq. To be specific, here we consider Mq = 10~ 5 vpq, 
and depict the resulting relic axino density in the left panel of Fig. for Tr > 10 3 GeV 
and vpq = 10 11 GeV (solid), 10 12 GeV (dashed), and 10 13 GeV (dotted). For numerical 
result, we use g 2 = g 2 (T = 10 6 - 10 8 GeV) ~ 1, and find 



Ya{T R > Mq) ~ 8.8 x 10 



-d 



10 12 GeV\ 2 / 


' Mq * 


vpq J \ 


,10 7 GeV, 


10 12 GeV\ 2 


( T R 


VPQ J 


V10 7 GeV 



y a (T<0.1M e) = 2.2xlO-(^) (3.16) 

Again, the black dashed line in the figure represents the axino number density when axinos 
were in equilibrium with the thermal bath of gauge-charged particles. Our result shows 
that axinos could indeed be in thermal equilibrium if vpq < 10 12 GeV and the reheat 
temperature is high enough. 

In the right panel of Fig. 1C, we plot the contours giving £lah 2 — 0.1, which shows that 
the CDM constraint gives a severe upper bound on the reheat temperature depending upon 
the axino mass. Note that our analysis includes the axino productions by the exotic quark 
multiplet Q,Q C , which in fact provide dominant production channels for T > O.IMq. As a 
result, for certain range of the axino mass, the upper bound on Tp can be more stringent 

than the bound one would obtain based on the effective interaction ^-4 f d 2 6AW a W a 

alone. Still, one can relieve the bound by assuming that Mq is small enough, e.g. close to 
the weak scale, for which the situation becomes similar to the case of DFSZ model. 

As in DFSZ model, we should consider the BBN constraint on the decays of long-lived 
NLSP. For the neutralino decays x ~^ l/Z + the rate is given by ( 3.14|) as in DFSZ 



model. However, there is no tree level coupling for the stau decay f — > t + a in KSVZ 
model, and the decay rate is suppressed by the two-loop factor as (L c 



r(r r + a) ~ — _ — — * In ' 



167T \87r 2 cos 4 6w J V rn 2 J \2v 2 P q 
1 / rw, \ 3 /l0 12 GeV x2 



10 4 sV200 GeVy V V PQ J 



(3.17) 



where a em is the fine structure constant and eQ = 1/3 is the C/(l) em charge of Q. Hence, 
neutralino lifetime is the same as in the DFSZ case, while stau lifetime is 0(1O~ 2 ) - O(10 4 ) 
s for vpq = 10 9 — 10 12 GeV. Therefore, the neutralino NLSP is safe as before, while the 
stau NLSP is safe only for a relatively heavy m? ■ 

4. Conclusion 

In this paper, we have discussed certain features of the effective interactions of axion su- 
permultiplet, which might be important for the cosmology of supersymmetric axion model, 
and examined the implication to the thermal production of axinos in the early Universe. 
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If the model has a UV completion at a fundamental scale M* 3> vpq, in which the PQ 
symmetry is linearly realized in the standard manner and all non-renormalizable interac- 
tions are suppressed by the powers of 1/M*, the axion supermultiplet is decoupled from 
the gauge-charged fields in the limit M^/vpq — > and vpq/M* — > 0, where M$ is the 
mass of the heaviest PQ-charged and gauge-charged matter multiplet in the model. As 
a result, in models with small values of M^jvpQ and vpq/M*, the axino production rate 
at temperature T S> M$ is suppressed by the powers of small M$/T. Such decoupling 
feature is not manifest in generic form of the effective lagrangian of axion supermultiplet, 
however it should be imposed as a matching condition at the scale just below vpq. 

Our observation is particularly important for the cosmology of supersymmetric DFSZ 
axion model in which M$ corresponds to the MSSM Higgs /i-parameter which is far below 
the PQ scale vpq. Cosmology of KSVZ axion model can be significantly altered also, if 
the PQ-charged exotic quark has a mass well below vpq. We have performed an explicit 
analysis of the thermal production of axinos for both the DFSZ and KSVZ axion models, 
and presented the resulting relic axino density as well as the bound on the reheat tempera- 
ture Tp coming from the observed dark matter density in our Universe. Our analysis does 
not take into account the effects of SUSY breaking and electroweak symmetry breaking. 
For Tp close to the weak scale, a more careful analysis is required, including the NLSP 
decays into axino as well as the mixing due to SUSY breaking and electroweak symmetry 
breaking, and it will be the subject of future work [31|. 
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Appendix A: 1PI axino-gluino-gluon amplitude 



In this appendix, we provide an explicit computation of the 1PI axino-gluino-gluon ampli- 
tude for a simple, but still general enough, form of Wilsonian effective interaction of the 
axion superfield at a generic cutoff scale A < vpq, which is given by 



AiC 
A 2 £ 
A 3 £ 



d 2 e Cw A 



32vr 2 



V PQ 



W aa W* + h.c, 



$Q + ^) ( J, ,., L .vviV'W J_ r, ,,T>t,T, , _i_ ,V<TV'"'<V 



d 2 e 



VPQ 

A 

VPQ L 



(xi + xf )Afi*i$; + (x 2 + x c 2 )M 2 ^2 



+ h.c, 



(4.1) 
(4.2) 
(4.3) 
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(d) (e) 
Figure 11: Axino-gluon-gluino amplitudes from the general effective interactions 



where W° are the gluon superfields, $i + <l>£ (i = 1,2) form 3 + 3 of SU(3) C , and M 2 > Mi, 
but both masses are well below the cutoff scale A. In this effective theory, the PQ symmetry 
is realized as 

A^A + ivpqa, $; -> e* 2 ****, -)• e i5 ? a $?, (4.4) 
and the corresponding PQ anomaly is given by 

Cpq = CV + J^(x< + x c i). (4.5) 

Under the field redefinition 

$i -»• e^/^o^i, $| -> e^/^«$f, (4.6) 
the Wilsonian couplings change as 

CV ^ C w + ^{zi + z c i) } 

i 

( yf , %i, %i) ^ (m + Zi, Vi + zt, Xi - Zi, x\ - zf). (4.7) 

It is straightforward to compute the 1PI axino-gluino-gluon amplitude from the above 
effective interactions, yielding (see Fig. 11) 

q 2 

A 1PI (k,q,p) = - y 5 A {k + q+p)C lPI {k,q,p)u{k)a llv ^v{q)e^p u , (4.8) 

16tt z \/2vpq 

where k, q and p are the 4-momenta of the axino, gluino and gluon, respectively, and 

cipi = c w - + Hi) + + $ + ^ + *f) F & v M i) ( 4 - 9 ) 
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for 



F(jp,q;M) = 



dx 



l-x 



dy 



2M Z 



M 2 — [p 2 x(l — x) + 9 2 y(l — y) + 2(p • <7)xy] 



(4.10) 



Note that C\pj is invariant under the reparametrization ( |4.7| ) as it should be. For the 
axino production by gluon supermultiplet, the relevant kinematic situation is that axino 
and gluino (or gluon) are on mass-shell, i.e. k 2 = q 2 = in the limit to ignore SUSY 
breaking effects, while gluon (or gluino) is off the mass-shell with p 2 = (k + q) 2 < 0. In 
such kinematic region, F(p, q, M) is given by 



F(k 2 = q 2 = 0,p 2 ^0) = ^ 

\P I 



log 



1 + v 7 ! -4M 2 /V| 
-1 + a/1 - 4M 2 /|p 2 



(4.11) 



which can be approximated as 



F 



M 2 



\p 2 



\p 2 



1-2,1^ T-o, for \p z \^M 2 



F ~ 1 



12 M 2 



for \p 2 \ < M 2 . 



We then have 
Cip/(A: 2 = 

with 



0, p 2 + 0) = C 1P7 + o 



p 



M heavy, 



light , 2 



(4.12) 
(4.13) 

(4.14) 



C 1P i(M$ < p 2 < A 2 ) = C w - (yi + y\ + m + ^f), 
C\pi(M\ < p 2 < M 2 ) = C W + (x 2 + x c 2 -yi- y c x ), 
Cipi{p 2 < Ml) = C w + (xi + + x 2 + 



(4.15) 



If &2 + ^2 corresponds to the heaviest PQ-charged and gauge-charged matter superfield in 
the underlying model, we would have the boundary condition 



C 1PI (M 2 <p 2 < A 2 



0, 



(4.16) 



which yields the following expression of C\pi at generic momentum scale below vpq: 



Cipi(p) = Yj & + S i +& + Vi 

M?>\pi\ 



(4.17) 



Finally, we consider the expression of C\pj in another kinematic situation when one 
of the 4-momenta p, q, k is vanishing: 



F(k = 0,p 2 = q 2 =£ 0) 



2M 2 /|p 2 



y/l + AM 2 /\ P 2 



In 



1 + v / l + 4M 2 /\p 2 \ 
-1 + a/1 + 4M 2 /V| 



(4.18) 
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which has the limiting behavior 



F 
F 



2M 



,, ln (^J fOT \P \» M ' 
for \p z \ < M l 



\v*\ \\v 



6M 2 



and therefore 



M heavy, 



P 



(4.19) 



, (4.20) 



where Cf PI are given as ( 4.1 5| ) . 



Appendix B: RG running and threshold corrections 

In this appendix, we discuss the RG running of the Wilsonian and 1PI amplitudes of the 
axion superfield, using the method of analytic continuation into N = 1 superspace. For 
this, let us consider a generic Wilsonian effective lagrangian 



C W (A) = J d 4 9Z n (A;A + A^l$ n 

+ / d 2 6 Iff (A;A)W aa WZ + ^e-^ + ^ A ^M mn ^ n 



1 

6« 



+ 7(; -(i m +J„+i P M/, P(3Amn!)$m$n$p + h c 



for which Cw and y n can be defined as 

Off 



C W {A) = -8tt z v pq 



dA 



A=0 



dlnZ n 

y n (A) = v PQ 



OA 



(4.21) 



(4.22) 



(4.23) 



A=0 



We then have 



dx n 
din A 
dCy/ 



0, 



" 87r v pq in 



d fdf 



ff 



din A dA VdlnA 

while y n can have a nontrivial RG running as 



A=0 



dy n 
din A 



d ( d In 2L 



<9A V din A 



(4.24) 



(4.25) 



A=0 



where we have used that the RG running of the holomorphic gauge kinetic function ff is 
exhausted at one-loop, and therefore dff/dlnA = 6 a /167r 2 is an A-independent constant. 

In the appendix A, we have seen that the 1PI axino-gaugino-gauge boson amplitude 
C^pj in the kinematic regime with k = and M 2 < p 2 < Mf is given by 



p 2 \ _ (Ml ( p 2 



c 1PI = c 1PI + o{ * + f 'ta I* 



(4.26) 
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where C\pi is a constant in 1-loop approximation. Including higher loops, C\pi can have a 
logarithmic p-dependence, which can be determined by the 1PI RG equation. To examine 
this, let us consider the 1PI effective lagrangian at a momentum scale p < A, including 



-ipi 



d 4 9 



Z n ( P , A + A*)$t $ n + — J" a (p, A + 



D 2 

W aa —W? + h.c. 

p2 



,(4.27) 



where Z n is the 1PI wavefunction coefficient chosen to satisfy the matching condition 

Z n (p 2 = A 2 ) = Z n (A), (4.28) 
where Z n (A) is the Wilsonian wavefunction coefficient at the cutoff scale A. We then have 



C\pi = -8tt 2 vpq 



dA 



A=0 



One can introduce also the canonical 1PI Yukawa couplings 

e -(Xm+Xn+Xp)A/v P Q ^ 



Amnp — 



^mnp 



which obeys 



vpq 



mnp\ 



OA 



A=0 



where 



\/ Z m Z n Z p 

= -(x m + x n + x p + y m (p) + y n (p) + y p (p))\A mnp \ 2 , 
d In Z n (p) 



Vn{p) = V P Q 



dA 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



A=0 



The 1PI gauge kinetic function obeys the Novikov-Shifman-Vainshtein-Zakhatov (NSVZ) 
RG equation 



din p 2 



where 



7n(p) = 



ciln Z n 
dlnp 



(4.34) 



Both /3 a and j n can be expressed as a function of the 1PI gauge coupling g 2 = l/J- a and 
the 1PI Yukawa coupling |A mnp | 2 . We then find the 1PI RG equations for C^ PI and y n are 
given by 



dCipi 
dlnp 



djjn 
dlnp 



a 2 d /3 a 
-Svr VPQ ^ 
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(4.35) 
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(4.36) 
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In section 2, we noticed that the Yukawa couplings of the model are constrained by 
the PQ selection rule 



y m + y n + y P )A 



mnp 



J PQ , 



(4.37) 



With the boundary value C^ PI at p ~ upq, which is estimated as ( p. 73 ), the RG equations 
(4.35) and ( 4.36| ) assure that 



C{ PI {M^ <p< v PQ ) = O 



Ml 



J PQ , 



+ o 



U PQ 
AG 



(4.38) 



and therefore the result (|2.66|) remains valid even after higher loop effects are included. 



Appendix C: Thermal axino production by matter multiplet 

In this appendix, we present a calculation of the relic axino density when the axinos are 
produced dominantly by the matter multiplet $ at T > M$. The Boltzmann equation of 
axino production is given by 

din ~ 

— + 3Hria = {v(i+j^K+a)v)(ninj - n K n~ a ), 



(4.39) 



where nj is the number density of the I-th particle species. At high temperature of our 
interest, all particles except axinos are in thermal equilibrium, so we can write ni = n eq , 
where is the equilibrium number density of particle /. From the relation of the detailed 
balance, we have 

{<J(l+j^K+~a)v)nf nj = {cr {I+ j^ K+ ~ a) v)n£nl q . 
Defining Y& = ria/s, the Boltzmann equation becomes 



^ = s(a {I+J _ K+ - a) v)Y^-Y a ). 



(4.40) 
(4.41) 



For A„ = Y~ — Ya, we have 



-s((T( I+ j^ K+ a)v)YK q dt, 



which yields 



A 5 (T ) = A a (T i? )exp 



Tr 



dT 



(4.42) 



(4.43) 



'T H(T)T 

As long as axinos are relativistic in the temperature of interest, Y? q is independent of 
temperature. Then we can use 
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(4.44) 
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where the relation ( 4.4CQ is used for the last identity. Under the initial condition Aq(Tr) 
Y- eq — Ya(Tji) = Y? q , the above result gives 



Y a (T ) = Y? 1 - exp 



1 [ Tr (^(i+J^K+a^nfnJ 1 



Tr 



To 

eq eq 



s(T)H{T)T 



dT 
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16vr 4 



dt t z K x {t) 
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tT R 
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(s — mj — rrij) — 4mjrrij 



where g = 135\/T0/(2vr 3 ^ /2 ), M Pi = 2.4 x 10 18 GeV, t R = (m/ + mj)/T R , and K^t) is 
the modified Bessel function |32[|. The cross section cr(s) can be obtained by 



da 



1 



, cift 7 2£ pi 2£ p2 \v pi - v P2 1 (2vr) 2 4£' c ___ 
with the amplitudes listed in Table [|. Then, for Tr S> M$ 3> M^, we find 
vAsCfrr, gM Pl f 8g 2 (x* +x^Ml\T? J \ 2 \ 1 



(4.45) 
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V 2vr 4 / V 



iV^-l 



(4.46) 



where \T^\ = Tr(T a T a ) = (N — l)/2 for the SU(N) gauge group, and the numbers in the 
first line come from the integrations for the processes C, D, E, G, I, H and J, respectively, 
in Table |l[ From the above result, the relic axino mass density can be determined as 



Qah 2 = Pah 2 /p c = man- a h 2 /p c = maYas(T )h 2 1 p c 



where we have used p c /[s{T Q )h 2 } = 3.6 x 10~ 9 GeV. 
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